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Abstract 

We present a new proof of Witten's conjecture. The proof is based on the analysis of 
the relationship between intersection indices on moduli spaces of complex curves and 
Hurwitz numbers enumerating ramified coverings of the 2-sphere. 



< 

1 Introduction 

Let M.g- n denote the Knudsen-Deligne-Mumford moduli space of genus g stable complex 
curves with n marked points For each i € {1, . . . , n}, consider the line bundle Li 
over A4 g - n whose fiber over a point (C;x\, . . . ,x n ) € -Mg;n is the cotangent line to the 
curve C at the marked point Xj. Let tpi £ H 2 {M. g - n ) denote the first Chern class of this 
line bundle, tpi = c i(A)- Consider the generating function 

F(to, h,...)= ^{r dl . . . r dn ) l%'- tdn , M (1) 
^ |Aut(di,...,d n )| 

\ for the intersection numbers of these classes, 

a 

(r dl ...r dn )= [_ 
' J Mg;n 

•i-H . 

^ , where the genus g is uniquely determined from the identity 

— 

d\ H h d n = dim .M g;n = 3g - 3 + n. 

The first few terms of this function are 

1 1,1,1 1 o 1 1*1 1 9 

24 6 48 1 24 6 1152 72 1 12 48 

1 3 2 1 4 29 1 1 

+ 6 tok + 24 ^ 2 + 5760^ 3 + 384 ht4 + llM^ 5 + ' ' ' 
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The celebrated Witten conjecture asserts that 

the second derivative U = d 2 F/dtf i of the generating function F satisfies the KdV equa- 
tion 

dU_ _ dU_ J_d*U_ 

dh ~ dt^ + l2~d^' ( ' 

The motivation for this conjecture can be found in |23| . and for a detailed exposition 
suitable for a mathematically-minded reader see, e.g., ^3]. Note that the KdV equation can 
be interpreted as a reccurence formula allowing one to calculate all the intersection indices 
provided "initial conditions" are given. Witten has shown that the function F satisifes the 
so-called string and dilaton equations, which together with the KdV equation generate the 
whole KdV hierarchy and provide necessary initial conditions. 

A number of proofs of the Witten conjecture are known, but all of them exploit tech- 
niques that do not seem to be intrinsically related to the initial problem: Kontsevich's 
proof makes use of Jenkins-Strebel differentials and matrix integrals, the proof due 
to Okounkov and Pandharipande |18| . which starts with Hurwitz numbers, also involves 
matrix integrals and graphs on surfaces study, as well as asymptotic analysis, and, finally, 
Mirzakhani's proof |15j is based on the Riemannian geometry properties of moduli spaces. 
The goal of this paper is to present a new proof using purely algebro-geometric techniques. 
Similarly to the proof due to Okounkov and Pandharipande, we start with Hurwitz numbers, 
but then we follow a different line. 

Hurwitz numbers enumerate ramified coverings of the 2-sphere with prescribed ramifi- 
cation points and ramification types over these points. We deal only with ramified coverings 
whose ramification type is simple over each ramification point but one. Our proof is based 
on the following, now well-known, properties of these numbers: 

• the ELSV formula [HIE], see a ls° [Hji |13j . relating Hurwitz numbers to the intersection 
theory on moduli spaces; 

• the relationship between Hurwitz numbers and integrable hierarchies conjectured by 
Pandharipande ^H] and proved, in a stronger form, by Okounkov 17 . 

Using the ELSV formula we express the intersection indices of the ^-classes in terms of 
Hurwitz numbers. The partial differential equations governing the generating series for 
Hurwitz numbers then lead to the KdV equation for the intersection indices. Note that the 
existence of such a proof has been predicted in |7j . One of the main features of the proof 
consists in the fact that known effective algorithms for computing the Hurwitz numbers, 
which are relatively simple combinatorial objects, lead to an independent tool for computing 
the intersection indices. We describe the properties of the Hurwitz numbers in detail and 
deduce Witten's conjecture from them in Sec. 2. Section 3 is devoted to a discussion of the 
proof. 

We express our gratitude to the participants of our seminar at the Independent Uni- 
versity of Moscow for useful discussions. Our personal gratitude is due to S. Shadrin and 
D. Zvonkine, whose papers contain, in an implicit form, the idea of inverting the ELSV 
formula. Special thanks are due to D. Zvonkine for careful proofreading of the first versions 
of the paper. The final version of this paper has been written at the Max-Planck Institut 
fur Mathematik, Bonn (MPI preprint 05-55), to whom we are grateful for hospitality. 
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2 Proof 



2.1 Hurwitz numbers 



Fix a sequence b± 



. . . , b n of positive integers. Consider ramified coverings of the sphere S 2 



by compact oriented two-dimensional surfaces of genus g with ramification type (b\, . . . ,b n ) 
over one point, and the simplest possible ramification type (2, 1, 1, . . . , 1) over all other 
points of ramification. According to the Riemann-Hurwitz formula, the total number m of 
these points of simple ramification is 



where B = b\ + • • • + b n is the degree of the covering. If we fix the ramification points in 
the target sphere, then the number of topologically distinct ramified coverings of this type 
becomes finite, and we denote by hg.^ ^ the number of these coverings, with marked 
preimages of the point of degenerate ramification, counted with the weight inverse to the 
order of the automorphism group of the covering. These numbers are called Hurwitz num- 
bers. 

The ELSV formula [U |H] expresses the Hurwitz numbers in terms of Hodge integrals 
over the moduli spaces of stable complex curves: 



for g > 0, n > 1 or g = 0, n > 3. The numerator of the integrand is the total Chern class 
of the vector bundle over A4 g - n dual to the Hodge bundle (whose fiber is the g-dimensional 
vector space of holomorphic differentials on the curve) , Xi G H 2i {M g - n ). The integral in © 
is understood as the result of expanding the fraction as a power series, with further selection 
of monomials whose degree coincides with the dimension of the base (there are finitely many 
of them) and integration of each of these monomials. 

The integral (jl} is a sum of intersection indices of both ip- and A-classes. In |18| . the A- 
classes are eliminated by considering asymptotics of integrals of this kind. In contrast, in the 
present paper, the exclusion of the A-classes is based on simple combinatorial considerations 
originating in [2HH1], see Sec. 2.2 below. 

Now consider the following exponential generating function for the Hurwitz numbers: 



where the summation is taken over all finite sequences 61,..., b n of positive integers and 
all nonnegative values of g, with m given by Eq. Q- According to Okounkov [17] . the 
exponent e H of this generating function is a r-function for the KP-hierarchy. In fact, 
Okounkov proved a much stronger theorem stating that the generating function for double 
Hurwitz numbers (those having degenerate ramification over two points rather than one) 
satisfies the Toda lattice equations. We do not need this statement in such generality, and 
in Sec. 3 we discuss a simple proof of the fact we really need. (Various kinds of Hurwitz 
numbers have been since long known to lead to solutions of integrable hierarchies, but we 



m = 2g - 2 + n + B 



(3) 



hgM,...,b n = m 




(4) 
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were unable to trace exact statements and origins.) The function e H being a r-function for 
the KP hierarchy means, in particular, that H satisfies the KP- equation 



d 2 H d 2 H 1 /d 2 H\ 2 1 d 4 H 



dp?, dpsdpi 2 \ dp\ ) 12 dpf ' 
We use this equation below to deduce from it the KdV equation for the function U. 



(5) 



2.2 Expressing intersection indices of ^-classes via Hurwitz numbers 

Obviously, for each nonegative integer d there exist constants cf, b = 1, . . . , d+ 1, such that 



EtZ^ = ^ + °(^ +1 )' (6) 

6=1 V 

and these constants are uniquely determined by this requirement. They are given by the 
formula 



c b 



(d-6 + l)!(6-l)!" 

Indeed, we need to prove that the first d — 1 derivatives in ip of the linear combination 



,, , . ¥> 

vanish at 0, while the d th derivative is d\. The i th derivative of this linear combination 
evaluated at ip = is 

The expression in brackets coincides with the result of applying the i th iteration of the 
operator xd/dx to the polynomial (1 — x) d and evaluating at x = 1, which is for < i < d 
and ( — l) d d\ for i = d. 

Multiplying identities © for different d we obtain the following equality: 

di + l d n +l di d n n 

V V bl ' ' ' bn = TT^ di + 



where dots on the right-hand side denote cohomology classes of degree greater than d\ + 
h d n . This means, in particular, that for d\ + • — h d n = 3g — 3 + n the linear combination 

di + l dn+l ,. 1 \ i J_ \ 

V r dl „d n / 1 ~ Ai H ± A g 

b1 "' bnJ ^ (i-w-a-w 

is simply (r^ . . . T^ n ), because the integral of the terms of higher degree vanishes. Taking 
into account the coefficient of the integral in Eq. (JIJ), we obtain the following explicit identity. 
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Theorem 2.1 For any sequence of non-negative integers d\, . . . , d n we have 



di+l d n +l ( 1 n , 1 N d . +1 _ 6 \ 



6 1= 1 6„=1 \ i=l V"" 1 1 / 

where g is determined by the left-hand side, dj = 3g — 3 + n, and m = 2g — 2 + B + n. 

It is convenient to reformulate the statement of the theorem in terms of generating 
functions. Decompose the generating function H into the sum 

H = Hq-i + Hq-2 + H st , (7) 

where the stable part H s t contains all the monomials whose coefficients are given by the 
ELSV formula Q, and H^i and -ffo;2 ar e the generating functions for the numbers of 
ramified coverings of the sphere by the sphere with 1 ("polynomial") and 2 ("trigonometric 
polynomial") preimages over the distinguished ramification point, respectively. The latter 
generating functions are known since Hurwitz E3 ; see also : 

00 gb-1 00 fob-2 
6=1 V ; ' 6=1 

1 ^ ^61+62 1 °^ b h b b2 

Ho;2 = j 22 h 0M2 P blPb = - >J + \ 2 " bllb2] Pb,Pb 2 P bl+b2 

61,63=1 ^ ' 61,62 = 1 V 

(note that this case can also be considered as been covered by the ELSV formula, but with 
the moduli spaces A4 g - n replaced by the moduli stacks A^q ; i, M.o-2)- In fact, we are going to 
use below not the precise formulas for iio;i> -£^0;2, but the fact that they contain only terms 
of degree at most 2 in p^, which yields 

d 2 1 
— (H 0; i + H 0;2 ) = -p 2 , 

d 2 

Tr-O (H0;1 + Hq-2) = (3, 

dpi 

ih -(H , 1+ H , 2 ) = V. 



dpidp 3 

Denote by G st = G st ((3; to, t\, . . . ) the result of the following change of variables in the series 

»- E ^^ r^t, (8, 

d=b-l y ' 

The result of this substitution is a series in to,t\, . . . whose coefficients are formal Laurent 
expansions m /3 2 / 3 . Indeed, the powers of (5 in the contribution of a monomial p bl . . . p bn to 
the expansion have the form 

n ( 2d- + 1\ 2 / n \ 

(61 + • • • + b n + 2g - 2 + n) + {~ b i J = 3 I % ~ 3 + n ~ Yl di J ' 

hence become even integers when multiplied by 3. On the other hand, the powers of [3 at 
each td are bounded from below, because ta enters only the expansions for p\, . . . ,Pd+i- 



5 



Theorem 2.2 1. The series G s t contains no terms with negative powers of p. 

2. The free term in f3, G st |/3 = [which is correctly defined due to the first statement) 
coincides with the generating function F for the intersection numbers given by Eq. (|T|). 

F(to, ti, . . . ) = G s t(0; to, h, . . . ). 

Proof. Collect together the terms of the series H s t corresponding to given values 
n, b\, . . . , b n , for all g, and set 

H bl ,...,b n = Yl h g]bu ..., bn — . 

m=B+2g-2+n 

Then the ELSV formula Q) can be conveniently rewritten as 

H _ .B+n/3 u b l / lz /^ 2/3 Ai + /? 4/3 A 2 - P^Xs + ■ ■ ■ \ 

bu-..,K P ii w y (1 _ ^2/3^) . . . (1 _ 6^2/3^) J > 

where we understand the numerator as the formal sum and the angle brackets mean inte- 
gration of each monomial over the space M g - n whose dimension coincides with the degree of 
the monomial. Indeed, consider a summand containing the integral of ip^ 1 . . . ifr^Xj in the 
expansion on the right-hand side. The genus g corresponding to the domain of integration 
Aig- n is computed from the equality 

^ dj + j = 3g - 3 + n. 

This summand contributes to a term of degree m in f3 iff the relation 

„ n 2 J—*- 2 „ 
m = B + - + - + -j = 2g - 2 + n + B, 

i=l 

which is exactly the relation between g and m in the definition of the series H bl holds. 

Now, the explicit form of the change of variables (JHJ) implies that the coefficient Gr^ ... <2„ 
of the monomial t^ ■ ■ ■ td„/\Aut(di, . . . , d n )\ in G s t is equal to 

di+l d n +i / n ^ k . x 2d i+ 1 \ 

= ■ • • n (j.-i.+ivft 4 *- 1 ^.-^ 

6l= i 6n= i Vi=i ( 1 l+ J 1 / 

SEif ^ ^ / eg . . . cfc (1 - Al + A2 _ ^6/3^ + _ ) V 

(l-&l/3 2/ ¥l)---(l-6n/3 2/3 ^n) / 

n 

= /H E * (^2/3^ +...)(!_ /3 2/3 Ai + ^4/3 Aa _ )} 
i=l 

ra 



where dots in the last two lines denote terms of higher degree in /3. 
Theorem 12.21 is proved. 
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2.3 Reduction to the KdV equation 

By Theorem 12.21 G s t is a power series in /3 2//3 , to, t\, ■ ■ ■ whose coefficient G\p=o of f3° co- 
incides with the function F. By definition, G s t is the result of the substitution (JHJ) to the 
series H st = H — Ho-i — Hq-2- Equation © together with the explicit formulas for Hq-\ and 
Hq-2 imply that H st satisfies the equation 

8 2 H st _ d 2 H st 1 fd 2 H st \ 2 l p2 d 2 H st 1 d 4 H st 



dp 2 dpsdpi 2 \ dp 2 J 2 dp 2 12 dp\ 
The change of variables Q results in the following change of partial derivatives: 

dpi dt ' 

op 2 otx dt 

9 = Q pu/sd + iy3 d + 9 10/3 d 

Op3 Ot2 Oti 2 Oto 

After substituting this into the equation above and dividing the result by /? 16//3 , we rewrite 
it as 

/ d 2 G st l(d 2 G«\ 2 1 <9 4 G st \ . „, /3 / a 2 G st d 2 G s 



\dtidt 2 V dt 2 J 12 dt% 
The coefficient of (3° in Q has the form 

d 2 F 1 fd 2 F\ 2 1 d 4 F 



0. 



Differentiating the last equation over to we obtain the desired KdV equation © for the 
function U = d 2 F/d 2 t 2 ) . This completes the proof of Witten's conjecture. 



3 Odds and ends 

Trying to make the present paper more self-contained, we discuss in this section the Hirota 
bilinear equations and the KP hierarchy, as well as an explicit presentation of the function 
e H as a r-function for the KP hierarchy. Although specialists in integrable hierarchies are 
well aware of these facts, it is not an easy task to find their compact and readable exposition. 
We refer the reader to [201, ffS]' 0> an d f° r a description of the relationship between 
integrable hierarchies and the geometry of semi-infinite Grassmannian. The KP hierarchy is 
a system of partial differential equations for an unknown function H. The KP equation © is 
the first equation in this system. Similarly to the case of the KdV equation, the expansion of 
a solution to the KP equation can be reconstructed from "initial conditions" . The exponent 
r = e H of a solution H to the KP hierarchy is called a r-function of the hierarchy. The 
equations of the KP hierarchy rewritten for r- functions also are partial differential equations; 
they are called the Hirota equations. They possess a nice property of being quadratic with 
respect to r. 
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3.1 Semi-infinite Grassmannian, Hirota— Pliicker 
bilinear equations, and integrable hierarchies 

Define the charge zero Fock space T as the completion of the infinite dimensional coordinate 
vector space over C whose basic elements s\ are labeled by partitions, 

Recall that a "partition is a nonincreasing sequence of integers A = (Ai, A2, . . . ), Ai > A2 > 
• • • > 0, having finitely many nonzero terms. Elements of T are infinite formal linear 
combinations of the vectors s\. We shall use the following two realizations of the Fock 
space. 

(1) The space J- can be identified with the space T = C[[pi,p2 5 • • •]] of formal power 
series in infinitely many variables pi,P2, ■ ■ ■ by setting s\ to be the corresponding Schur 
function. The Schur function corresponding to a one-part partition is defined by the ex- 
pansion 

so + siz + s 2 z 2 + s 3 z 3 + s 4 z 4 + ■■■ = e Pi*+P2#+P3^+... 5 
and for a general partition A it is given by the determinant 

s x = det ||sA 3 -j+i||- (10) 

The indices i,j here run over the set {1, 2, . . . , n} for n large enough, and since \ = for i 
sufficiently large, the determinant, whence s\, is independent of n. Here are a few first 
Schur polynomials: 

so = 1, s 1= pi, s 2 = ^(pl+ P2 ) i s 3 = ^(pl + 3p 1 p 2 + 2p 3 ), 

si,i = ^(Pi-P2), s 2 ,i = -{pf-ps), si,i,i = -{pi - 3p!p 2 + 2p 3 ). 

(2) Let V = C[z,z _1 ] be the ring of Laurent polynomials in z, which we treat as 
the vector space with the basis z', i £ Z. Identify T with the semi-infinite wedge space 
T = A°°/ 2 F freely spanned by the formal infinite wedge products of the form 

s x = z kl A z h ' 2 A z h[i A . . . , h=i- A,, 

for all partions A. Sequences fcj appearing on the right-hand side can be characterized as 
arbitrary strictly increasing sequences of integers satisfying ki = i for i large enough. 
The theory of the KP hierarchy can be summarized as follows. 

The projectivization PT = PA°°/ 2 V is the ambient space of the standard Pliicker 
embedding Gr PJ 7 , where Gr = Gr^^iV) is the Grassmannian of "half-infinite dimen- 
sional subspaces" , often referred to as the Sato Grassmannian. By definition, the elements 
of Gr are subspaces spanned by linearly independent vectors (p±, <f2, • • ■ in (the formal 
completion of) V such that for i large enough we have <pi = z l + . . . , where dots denote 
terms of lower order in z. Such a vector space can be interpreted as the wedge product 

r = ipi A (p 2 A f 3 A . . . (11) 
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Indeed, another choice of a basis does not affect this wedge product, up to a scalar factor. 
If all the functions <fi are Laurent polynomials and <pi is simply the monomial z l for all i 
sufficiently large, then the wedge product r can be represented, after expanding the brackets, 
as a finite linear combination t = ^ x c\s\, c\ £ C If the functions <pi contain infinitely 
many terms, then the function r is a formal linear combination of s\ and can be obtained 
in the following way: 

when expanding the brackets in the infinite wedge product pick one monomial sum- 
mand in each ifi in such a way that this summand is z l for all but finitely many indices i 
and do this in all possible ways. 

More explicitly, if <pi = Yljez a iJ z ^ t nen 

t = ^Jdet \\aij-Xi\\i,j>isx = detW}^ aj,kSj-k\\i,j>i- (12) 

Theorem 3.1 ( |20l 121 l^lj^ A (non-zero) function r is a r -function for the KP hierarchy 
iff the corresponding point [r] € PT belongs to the Grassmannian Gr C PT ' . 

In particular, each Schur polynomial s\ and any linear combination of the Schur func- 
tions Si corresponding to one-part partitions produces a solution to the KP equation (0). 

The image of the Pliicker embedding of the Grassmannian is known to be given by 
a system of quadratic equations. In our case, these algebraic equations for the Taylor 
coefficients of the function r take the form of partial differential equations for this function, 
known as the bilinear Hirota equations. 

3.2 Formulas for the generating function for Hurwitz numbers 

The exponent e H of the generating function for the Hurwitz numbers is nothing but the gen- 
erating function for the numbers of ramified coverings of the 2-sphere by all, not necessarily 
connected, compact oriented surfaces of Euler characteristic 2 — 2g. Take such a covering 
and let a point of simple ramification in the target sphere tend to the point of degenerate 
ramification. Then one can express the number of such coverings as a linear combination of 
the numbers of similar coverings with fewer points of simple ramification. This reccurence 
relation (the "cut-and-join equation" of [Hj) expressed in terms of generating functions reads 
as follows: 

Q e H i ™ / Q Q2 \ 

W = 2 £ [ {l + j)PiP ^ + ^ 'dp-dp- ) ^ = M ^ (13) 

where we denote by A the differential operator on the right-hand side. For an algebro- 
geometric interpretation of the cut-and-join equation, see \'22\ . 

Equation (|13|) can be solved explicitly. The operator A acts linearly on the space of 
weighted homogeneous polynomials in the variables pi,p2, ■ ■ ■ , with the weight of the vari- 
able pi equal to i. Moreover, it preserves the weighted degree of the polynomials, whence 
can be split into a direct sum of finite dimensional linear operators. The Schur functions s\ 
are eigenvectors of A, and they form a complete set of eigenvectors. 
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Denote by /(A) the eigenvalue of the eignevector s\. It can be checked that 

1 00 

/( A ) = 2E A ^- 2i + 1 )- 

i=l 

It follows that any solution of Eq. (|13j) can be represented as a sum 

A 

over all partitions A, for some coefficients c\. For the solution e , these coefficients can be 
computed from the initial value H\p = Q = pi, which yields 

e H = Y,sx(hO,0,...)s x e^. (14) 

A 

In ^7j, this form of the function e H was deduced from the representation theory of sym- 
metric groups. Note that c\ = (dim R\)/\\\\, where |A| = Ai + A2 + • . . and R\ is the 
irreducible representation of the symmetric group S\\\ corresponding to the partition A. 
Taking the logarithm we obtain a few first terms in the expansion of H: 

H = p 1 + I( e ^-2 + e - /3 )p? + -(e /3 -e- /3 )p 2 + — (e 3/3 -9e /3 + 16-9e- /3 + e- 3/3 )p? 
4 4 36 

+ ^ - 3e" + 3e-? - e" 3 > lP2 + ^ - 2 + e " 3 > 3 + . • • 
12 lo 

3.3 The exponent of the generating function for Hurwitz numbers as an 
element of the semi-infinite Grassmannian 

In Sec. 2.2, the fact that the function e H is a T-function for the KP hierarchy was established 
by a reference to Okounkov's paper ^7]. Here we present a more direct argument. 

For [3 = 0, the function e H = e Pl is the r-function corresponding to the semi-infinite 
space with the basis tfj(z) = z J e z . In other words, it is the result of the application of 
the operator e z 1 (the exponent being understood as the operator of multiplication by z^ 1 ) 
to the wedge product z A z 2 A z 3 A . . . (the vacuum vector) . 

Indeed, for this basis ipj, the coefficient of s\ is nothing but 5^(1,0,0, . . . ) because of 
the defining equation (Jlfljl and the relation Sfe(l, 0, 0, ... ) = l/k\ for the one-part partitions. 

The most natural approach to description of the evolution of this semi-infinite space 
along P consists in multiplying the coefficient of z l in each <pj by e'^ -1 ^' 2 . In other words, 
we apply to this space the exponent of the diagonal matrix with the diagonal elements 
i(i — 1)0/2. This procedure produces a well-defined basis in the /3-evoluting space. For 
this basis, however, Eq. (|12|) cannot be applied directly. To avoid this difficulty, one can 
normalize the basic vector tpj by dividing it by e- ? ^ -1 ^/ 2 thus making it start with 

^•( z ) = E(— -yr e 2 z - ( 15 ) 

Now it is obvious that Eq. ()12|) gives the coefficient of s>(l,0, 0, . . . )s\ equal to e-^ A ^, as 
desired, and we obtain 
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Theorem 3.2 The function e is the tau-function for the KP hierarchy corresponding to 
the semi- infinite subspace with the basis 1)15)) . 

This yields an independent proof of the following 

Corollary 3.3 The function H satisfies the KP hierarchy, in particular, it is a solution to 
the KP equation ©. 

The KP hierarchy degenerates into the KdV hierarchy if we are looking for solutions 
independent of variables with even indices (that is, the derivatives with respect to these 
variables vanish identically). In the normalization of the present paper, the variables tj are 
related to the variables pj by ti = j^zr\w_ P2i+i which, in particular, makes F independent 

of P2>P4) (Recall that (2i — 1)!! denotes the product of odd numbers from 1 to 2i — 1, 

(— 1)!! = 1.) Since the function F given by Eq. (^Q) satisfies the KdV hierarchy, its exponent 
e F is a r-function for the KdV hierarchy. Its explicit presentation similar to that of e H is 
given in [11| . It would be interesting to deduce the whole KdV-hierarchy for F from the 
KP-hierarchy for H directly, without referring to the string and the dilaton equations. 
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